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back-reaction is weak, evaporation proceeds via an adiabatic evolution, whereas for strong 
back-reaction, the decay proceeds in a somewhat surprising manner. Although information 
loss is inevitable in these models at the semi-classical level, it is rather benign, in that the 
information is stored in another asymptotic region. 
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1. Introduction 

One of the authors (M.O'L.) and T.Banks ^ had previously proposed that for a large 
class of modified scalar-gravity theories in which the classical geometries are all nonsingular, 
with causal structure identical to that of Reissner-Nordstrom, the Hawking evaporation 
to a final zero temperature remnant-like object could be studied without singularity, as 
opposed to the original CGHS (Callan-Giddings-Harvey-Strominger) models in which 
a now well known singularity was found. 

Here we report on calculations that support this picture. When in-falling matter 
perturbs one of these extremal solutions, two apparent horizons form. As the evaporation 
takes place, these apparent horizons approach each other. We find two distinct regimes, 
depending on whether the back-reaction is weak or strong in a certain sense. With weak 
back-reaction, an adiabatic approximation gives a correct description, and things settle 
back down to a stable remnant, with the apparent horizons meeting only after an infinite 
proper time. In the strong back-reaction regime, the apparent horizons meet after a 
finite proper time, and only after meeting do things settle back down to the extremal 
solution. Black holes in these models therefore evaporate in a completely nonsingular 
fashion, realizing the original objectives of CGHS. Information loss occurs at the semi- 
classical level, but only in a rather benign way. 

In section 2 we introduce the models of interest. We discuss in some detail the behavior 
near the double horizon of the extremal static semi-classical space-time in section 3. In 
section 4 we describe the adiabatic approximation for the nonsingular models. Section 
5 contains the results of our numerical analysis and section 6 is devoted to our conclusions 
and a discussion of their implications. 
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2. The Models 

Consider a Lagrangian taken from the general class of two-dimensional renormalizable 
generally covariant field theories 

= ^f^{D{4>)R + GmV4>f + . (2.1) 

We require that the potentials behave asymptotically like those of linear dilaton gravity 

0, 

as — > — oo. 

The renormalization group equations are hyperbolic on the two-dimensional target 
space of this model, and thus given a set of initial data one can consistently renormalize 
the model 0. In the following, without loss of generality, we will restrict attention to 
the class of models satisfying G{(j)) = ~2D'{(p). Other models may be obtained by a field 
redefinition of (p. Performing a Brans-Dicke transformation on the metric g = e~'^^g this 
Lagrangian may be rewritten in the simple form 

Li = ^g{D{^)R + Wm (2.3) 

where we have defined W{(j)) = e^'^H^cf)). This form of the Lagrangian is convenient for 
finding the classical solutions as described in [Q , in which reference the extended space- 
time geometries are discussed. 

All solutions are causally related to the two-dimensional r — t section of the four- 
dimensional Reissner-Nordstrom (RN) black hole (see fig. 1 for the extremal geometry), 
but as noted above, the geometries in these generalized models may be nonsingular. The 
nonsingularity is achieved by requiring that as (/> ^ oo, ~ e^^ and W ~ e^'^ , such that 



s 

D 




Fig. 1: Causal structure of the extremal zero temperature space-time ge- 
ometries. The arrows show the flow of time. CD is the Cauchy horizon for 
regions I and III. The global horizons are at fixed r and the geometry near 
r — Th is the same in both the DW and RN cases (in particular the distance 
to P in the direction of the arrow is infinite). The line ST is singular and at 
finite spacelike distance for RN, nonsingular at infinite spacelike distance for 
DW. 

n > m — 2. The multiple horizon structure is obtained by requiring that W have a zero. 
Henceforth this class of models will be referred to as DW models. 

Just as in four-dimensions one believes that if charge cannot be radiated away then 
the zero temperature extremal RN space-time with = Q^, where M is the mass and 
Q is the charge, will be an end-point for Hawking evaporation, it was conjectured that the 
zero temperature extremal limit of the DW models would also be a natural end-point for 
Hawking evaporation. 

To compute the back-reaction we look at the trace anomaly produced by the quantum 
fiuctuations of some conformally coupled (hence massless) matter fields. This is just the 
quantity calculated by Polyakov. We follow CGHS and study the back-reaction corrected 
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equations in the semi-classical limit, obtained by adding to the equations the terms corre- 
sponding to the trace anomaly of a large number of the matter fields fi and then taking 
to infinity with Nh fixed. Explicitly, the semi-classical action is 

N 

S = I (fx Cd-2Kd+pd-p + Y.^+f-^-f^ ' (2-4) 

in conformal gauge ( g±± = 0, g±^ = —^e^P, x^ = x*^ ± x^). Here we have defined 
K = Nn/12. 

By our choice of potentials, together with the condition D' + 2k < for all cp, we have 
avoided the singular kinetic term in field space which has been shown to be responsible for 
the singularity of the CGHS model. 

Let us note in passing that the general solution to the CGHS equations that is static, 
with mass above that of the vacuum, actually has a weak coupling singularity on the outer 
horizon (the quantum kink). This singular behavior arises through the interaction 
between the assumed dilaton-gravity asymptotics of the potentials (|2.2|) , and the large 
corrections. That is, this singularity is a weak coupling singularity that arises due to the 
assumed linear dilaton asymptotics of the Lagrangian, and so also will be present in our 
models. 

Note that we may alternatively require the potentials to behave asymptotically as the 
spherically symmetric reduction of the four-dimensional Einstein equations, 

D{(p) e-^-^, and H{(P) ^ 2 . (2.5) 

Here the two-dimensional metric is related to the four- dimensional metric via 

(")ds2 = ds2 + e-2W . (2.6) 

An example from this class of models, corresponding to the r — t section of the four- 
dimensional RN black hole will also be considered in the following, with the same one-loop 
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back-reaction term included as above. This model has recently been considered in 
Unlike the DW models described above, this model is singular at the classical level and 
has a large singularity at e~^'^ = k/2. 



3. Static Solutions of the Generalized Models 

The static solutions to the quantum equations of the DW models fall into three cat- 
egories, those that qualitatively have the same structure as the quantum kink of early 
studies of the CGHS Lagrangian, an extremal state that approaches the classical extremal 
state (as k — > 0), but with a singularity in the second derivative of the curvature on the 
horizon, and negative mass solutions that have the same structure as the classical negative 
mass solutions, in particular they are nonsingular. In the following we will refer to the 
extremal state as the vacuum of the two-dimensional theory, and measure mass relative 
to this vacuum state. 

Thus at this level we have found that quantum corrections have made the classical 
phase space more singular. However, the vacuum state is nonsingular and given that 
the positive mass solutions will radiate, it is plausible that none of the positive mass 
singularities will appear in the dynamical collapse and evaporation. 

In the next section we will introduce the adiabatic approximation used by Strominger 
and Trivedi to study the formally similar problem for a dimensional reduction of Reissner- 
Nordstrom geometry. The adiabatic approximation can be most easily investigated in 
light-cone gauge, so we now record the equations in that gauge. 

We have the line element. 



The scalar curvature is given by i? = —d'^h. To find the needed parts of the stress tensor. 



hdv^ + 2drdv . 
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we integrated the Bianchi identities V^Tq,^ = 0, using the equation for the trace anomaly 

i, 

T^"^ = 2T,^ + hr;!^ = '^R . (3.2) 

The components of interest to us are the trace, and l^'l^Tc^p^ where / = (/i/2, 1). 
The 4> and trace equations may be written as, 

2 _ {e-''^{D'W + k{W - \W')) - D"{D' + K){V(t)f) 



D'iD' + 2k) 



(3.3) 



2, {-e-^H'iW + W') + 2D"{S/ct^f) 
R = -9rh = ^7^^^^ , (3.4) 

where V'^cp = dr{2d^(p+hdr(j)), and {V(p)'^ = 2dr(pdy(p+h{dr(p)^ . The remaining local linear 
combination of the constraint equations is (conveniently the linear combination needed for 
matching across the shock wave), 

rfT^p = 2D\d^(t> + \hdr4>f + {\hdr + a, - \{drh)){\dr + d^)D 

K h 1 (3.5) 

+ ^{^dlh + - -^{drhf) + tr{r) 

where tr is determined by the initial state of the quantum vacuum, and for our collapse, 
and for static solutions with no net flux at infinity, we set tr to zero in coordinates that 
are asymptotic to the linear dilaton vacuum. 

For static solutions we put h = h{r), cj) = (j){r). With linear dilaton asymptotics at 
(j) — oo, W 4, D e"^"^, we find to lowest orderi, 

h{r) = 1 - 2Me-2'^ 

(3.6) 

(p{r) = — r . 

One can numerically integrate the static equations using the above as boundary conditions 
at r ^ oo. 



^ Certain choices of the potentials D and W will give an additional term in the asymptotics 
of h proportional to rexp(— 2r). 
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Now let us consider the behavior of fields near the horizon, r = r/^, of the extremal 
solution. We know from the classical analysis in that the classical extremal solution 
has to lowest order cj) = cpQ and h = hox^ near the horizon, where x = r — Vh- Let us look 
at the quantum corrections to these formulae, 

h = a\x^ + a2X^^'^ 

Plugging this into the static equations we find (evaluating all functions at ^ = ^o) 

W 

D'W = k{— - W) 



(3.7) 



(5 + l)(5 + 2) 



(X^ = 

T] = 5 

D"W + D'W - k{\W" - W) 



(3.8) 



D'{W + ^W') 
(3 and a2 are related by 

^ -Ka2{S + 2) 
^ 2ai6D' ^ ' 

Notice that the solutions have two obvious extensions through ^ = (see for a 

similar discussion). To study perturbations by shock waves of matter we will restrict to 

the odd extension which is the smoother of the two continuations and is the one that is a 

deformation of the classical geometry, 

(j) = (f)Q -\- [3x\x\^ 

(3.10) 

h = aix^ + a2x\x\'^^^ 

This also means that immediately above the shock wave we will find two apparent horizons 
at (f)j^ < (po and at (p- > (po, with a geometry qualitatively the same as the positive mass 
classical solutions. 



We can evaluate the parameters in the expansion near the horizon of the extremal 
solution and observe that the quantum vacuum near the horizon is indeed a small k de- 
formation of the classical vacuum. Let us fix 

D = e-^'^ -'y^e^'^, and W = 4-^'^e^'^. (3.11) 

For small k and large we find, 



2k (3 1 ■ 
ai = 4, d = — , — = - 
H a2 o 



(3.12) 



For definiteness we put (3 — —l/fJ. in the following. 

To check that the quantum behavior near (p = (po that we have displayed is consistent 
with the linear dilaton at infinity we numerically integrated out from (/> = to (/> = ±oo, 
and indeed have observed the linear dilaton vacuum for (p ^ —oo and the large cp classical 
behavior in the other asymptotic regime. 



4. The Adiabatic Approximation 

In order to study the semi-classical stability of the static extremal solutions of the 
DW models, consider sending in a matter shock wave 

/MPT/^ = 2Md{v). We define 
E = 2dy(p + hdr(p , (4.1) 

so that the future apparent horizon is the locus of E = 0. The discontinuity in E across 
the shock is 

= , . (4.2) 

To obtain explicit expressions we use the potentials of equation ( ^.11| ) and make a 
large small k, and small M expansion. Specifically, we have M/ fj, <^ k/h <^ 1. We may 



calculate the positions of the apparent horizons which turn out to be 



^.M, 1 , , 

r±=rh±y/JI{ — )^ (4.3) 
A* 



which is obtained by setting E + = 0. This implies that 



8 M 1+5 

a.E(r±) = T — ( — )^ . (4.4) 



The discontinuity in h across the shock is given by 

8M /• , , _ D' 



5h 



K 



dri-l-— ). (4.5) 

J ^ JD'{D' + 2k)' 



'D'{D' + 2k) 

By inserting this expression into the constraint equation we find 9„E at the apparent 
horizons, 

= ^ . 4.6 

One may then try to make an adiabatic approximation to compute the relative posi- 
tions of the horizons as a function of f , by assuming the equations (|4.4| ) and ( [4.6| ) continue 
to hold for all v if things are changing slowly enough: 

d.r^ = -^ = T-{r^-rn) (4.7) 

which gives 

f±(t;) =r;, + r±e-^("-"°) (4.8) 

or expressing things in terms of an effective mass, which measures the difference between 
the actual mass and the extremal mass 

M(^;) = Moe-i?^"-"°) . (4.9) 

The adiabatic approximation will break down when the energy fiux due to Hawking radi- 
ation is comparable to M{v), so will be valid as long as k/ /i <^ 1. 
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These expressions are almost identical to ones obtained in in the case of the spheri- 
cally symmetric reduction of the Reissner-Nordstrom black hole, apart from slight changes 
in numerical coefficients. The RN case is obtained by picking different potentials D and 
W . In this case, the adiabatic approximation is valid as long as n/Q"^ ^ 1, where Q is the 
charge of the extremal black hole. 



5. The Numerical Collapse 

When an extremal solution is perturbed by an incoming flux of matter the picture 
one expects is the following: above the shock wave two apparent horizons form on either 
side of the line and as the black hole evaporates these horizons approach each other, 
with the solution settling back down to the extremal solution at X'^. In this section we 
describe numerical solutions for a shock wave of in-falling matter impinging upon the 
extremal solution of the DW model described in section 3. We also consider the analogous 
calculation for the RN type model. Our aim is to test whether the extremal black hole is 
stable, and in what manner the evaporation proceeds. 

Another approach would be to do a linear stability analysis of the problem, where one 
might hope to be able to understand things analytically. In fact, it seems this is not so for 
the back-reaction corrected equations. A numerical analysis of this problem is tractable, 
but is of comparable difficulty to numerically solving the full nonlinear equations, so we 
choose to do the latter. 

In the context of the CGHS model, numerical results have previously been obtained in 
1^,11,12] . Here we are solving equations of precisely the same form, but with somewhat 



different potentials. See lO] for a description of the numerical algorithm used in this paper. 



and also appendices A and B, where the gauge and coordinate choices, and equations of 
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Fig. 2: Plot of motion of apparent tiorizons for DW model, in the strong 
baclc-reaction regime. The lower line is the outer apparent horizon, the 
middle line is (p = (po, and the upper line is the inner apparent horizon. 

motion are stated. For numerical purposes it is convenient to work on a grid of null 
lines, so conformal gauge is appropriate. This should be borne in mind when making 
quantitative comparisons with the results of the previous section, where it was necessary 
to use light-cone gauge to get analytical results. 

For the DW model of equation , with k = 10, = 15, 7 = 8, and shock mass 

Ms = 1.5, the results are plotted in fig. 2 . Note that here ^//i > 1 so one is in the 
strong back-reaction regime and the adiabatic approximation of the previous section is 
not expected to hold. The integration is stopped near x~ = where the line (p = (po starts 
out. This means these calculations will hold for either of the extensions described in the 
previous section. The horizons meet at finite proper time, and the line (pg goes from being 
spacelike to timelike at this point. Note the equations of motion imply that d^d-(p = 
when an apparent horizon ( where d+(p = ) intersects (p = (pQ. The position of the 
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Fig. 3: Plot of the curvature along a line of constant cj) — —1.2, for the DW 
model, in the strong back-reaction regime. Here the range of x'^ extends 
much further out. The curvature approaches a constant value which is the 
same as the initial curvature of the extremal solution. 



apparent horizons x {x'^) may be found by solving 

dx- _ dl(t> 



(5.1) 



dx+ d-d+(j) 

which indicates that dx~ /dx'^ blows up as d+d-cj) — > 0, as long as d^cj) remains finite. This 
is precisely what is happening at the meeting point of the apparent horizons, consistent 
with the numerical results. Also note that the wiggles in the path of the outer horizon 
are due to d-d+(f) becoming very small in this region. This causes the error in the path of 
the line d+(f) = to he much larger than, say, the error in the path of a line of constant 
(p. The conclusive indication that the horizons meet comes from the observation that the 
line (j) = (j)Q becomes timelike. Numerical convergence has been checked by varying the 
stepsize, and varying the position of the initial surface. 

Following a line of constant as x'^ becomes large, one finds the curvature approaches 
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Fig. 4: Plot of the path of a hne of constant (p — —1.2, for the DW model, 
in the strong back-reaction regime. As x'^ becomes large the line appears 
to asymptotically approach a null line. 

a constant which equals the initial curvature of the extremal solution, as shown in fig. 3 . 
This indicates that despite the fact that the apparent horizons have met, the solution still 
settles back down to the zero temperature extremal state. 

Lines of constant appear to approach a null line x~ = , which becomes a global 
horizon, as shown in fig. 4 . It is difficult to tell from the numerics whether Xq = 0, or 
whether it is shifted out to more negative x~ by the in-falling matter. More sophisticated 
numerical calculations are needed to answer this question definitively. 

Simulations were also performed for the weak back-reaction regime where k/ij, <S 1, 
when the adiabatic approximation is expected to hold. Here potential numerical errors 
were somewhat larger, but the results were found to be consistent with the adiabatic 
approximation. The apparent horizons persisted for very large x~^ , approaching the critical 
line = 00- 
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Fig. 5: Plot of motion of apparent liorizons for RN model, in the strong 
back-reaction regime. The lower line is the outer apparent horizon, the 
middle line is = 0O) and the upper line is the inner apparent horizon. 

Qualitatively similar results are obtained for the case of the RN model. In the strong 
back-reaction regime where k/Q^ > 1, with n = 200, = 60, and shock mass Mg = .08, 
results are shown in fig. 5 . One difference here is that coordinates are chosen so the line 
(j) = starts out at x~ = +oo. Again the solution appears to settle back down to the 
extremal solution as x'^ becomes large, after the apparent horizons have met. In the weak 
back-reaction regime {k/Q^ -C 1), results consistent with the adiabatic approximation 
were found, with the apparent horizons approaching the the critical line out to large values 
of x~^. 

6. Discussion and Conclusions 

The picture of the end-point of the evaporation of these two-dimensional extremal 
black holes which is consistent with the numerical calculations is shown in fig. 6 . In 
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Fig. 6: Penrose diagrams sliowing tlie end-points of tlie black liole evapora- 
tion. The case of strong back-reaction is shown on the left, the case of weak 
back-reaction is shown on the right. 



the strong back-reaction regime, lines of constant approach a global horizon, and the 
solution settles back down to the extremal one. After the apparent horizons meet the line 
= 00 turns timelike and becomes asymptotically null as it approaches the global horizon. 
In the weak back-reaction regime the apparent horizons meet at T^, and the critical line 
(f) = 4)0 remains spacelike. 

Our results indicate the presence of two qualitatively different regimes separated by 
some critical value of k/// in the DW case, and k/Q"^ in the RN case. This should not be too 
surprising, since a similar phenomena occurs in the case of the damped harmonic oscillator, 
where there is a critical value of the damping separating two regimes of behavior. Although 
we have studied two particular models, we believe that this behavior will be generic to the 
class of models in which W{(f)) has a simple zero, which includes the nonsingular models 
discussed in section 2. 
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Because for large k,/ n the apparent horizons coUide, locaUy, the end-point of the 
evaporation of these extremal solutions looks like an < static solution in the 
case of the RN model. One may wonder if some two-dimensional quantum positive mass 
theorem prevents such an unexpected occurrence. However, positive mass theorems only 
give us information about the asymptotic structure of the space-time and certainly do 
not preclude the type of relaxation that we observe here. Classically the DW and RN 



models obey positive mass theorems with the mass defined on space-like surfaces that 
become asymptotically null along and have their left boundary at the meeting of the two 
global horizons of the extremal solution below the shock wave. Unfortunately the quantum 
positive mass theorems do not appear to give any useful information in their present 
formulation due to the quantum corrections to the vacuum. Further, quantum positive 
mass theorems require a nonsingular spacelike surface that is asymptotically linear dilaton 
vacuum at both ends, see [|n|,|T5|,ff6l . In our situation this of course can never be the case, 



since the left boundary of our spacelike hypersurfaces meets the horizon of the extremal 
solution which looks distinctly not like the LDV. 

These results appear to confirm the conjecture that zero temperature extremal (in 
the sense of RN) two-dimensional black holes are semi-classically stable. In the case 
of the DW models the evaporation of nonextremal black holes proceeds in a completely 
nonsingular way, realizing the original objectives of CGHS. Information loss is inevitable 
in these models at the semi-classical levelS, since an in-falling flux of matter will always 
produce correlations with the region behind the global horizon, and hence inaccessible to 
an observer who finds themselves at ij^ after an infinite proper time. The information 
loss is of a somewhat benign type though. For observers outside the global horizon the 



See [|i|j3|,|l7| for related discussions. 
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quantum mechanics that they participate in is unitary, indeed the spacetime before the 
Cauchy horizon can (by definition) be fohated by space-hke Cauchy hyper surfaces. The 
information that is no longer accessible to them is stored in a stable remnant, one of 
an infinitely degenerate set of possible final states. This infinity of states corresponds to 
all the possible field configurations on a space-like slice through region III (fig. 1). These 
remnants avoid the problems of overproduction in external fields and divergences in virtual 
loops by nature of their large internal geometry . 

Let us note that the above discussion is based on semi-classical reasoning. It is still 
possible that when space-time is properly second quantized the information loss problem 
will be cured. Highly nonlocal quantum gravity effects may wind up giving the rem- 
nants a very long, but finite lifetime. If this lifetime is of order the age of the universe, 
experimenters making measurements over shorter times will still see an effective loss of 
information, which they would attribute to the existence of "stable" remnants. 

The conjectured extremal state possesses a Cauchy horizon as shown in fig. 1. We 
would like to comment on the sense in which this horizon is traversable!! and the possibility 
of information loss for observers who traverse it. We should first note that the Cauchy 
horizon in our model appears to be a double horizon, and as discussed in [§] has a softened 
divergence of the stress tensor as compared to the inner horizon of a nonextremal space- 
time. It is also possible that the disturbance produced by the shock wave will separate 
the local horizon and the global Cauchy horizon in a manner similar to the mass infiation 



models of Poisson and Israel \19 . In mass infiation it was shown that an infinite tidal force 



appears along the Cauchy horizon, but the singularity is weak in the sense that an 
observer can traverse it without getting stretched infinitely. 



See [18 1 and references therein for earlier discussions. 
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Let us assume then that a traversable Cauchy horizon is present. An observer who 
passes into the region above this Cauchy horizon faces a potential loss of unitarity asso- 
ciated with the lack of a global foliation of the space-time by space-like Cauchy hyper- 
surfaces. The existence of a Cauchy horizon, it's nature if it exists and the necessity of a 
unitarity restoring mechanism for those who cross it, are the subject of our ongoing inves- 
tigations. The answers to these questions should shed light on an S-matrix description of 
quantum gravity in the presence of remnant-like objects. 
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Appendix A. Numerical method for DW model 

The numerics are performed in conformal gauge 
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9++ = 9— = 0, 



9+- = 9-+ 




(A.l) 



The linear dilaton vacuum in the coordinates we choose is 



(j) = p = —^log{—x'^x ) 



(A.2) 
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while a solution with mass M looks like 



(j) = p = -^logiM -x+x~) (A.3) 
in the limit x~^x~ —oo. The equations to be solved are 
d+d-p -- 



d+d. 



D' + 2k 

D"d+(j)d-(j) + le'^P-^'^W + Kd+d-p 



D' 

Here D = e~'^^ — j'^e'^^ and W = A — /i^e^*^. The boundary conditions are that along 
x^ — 1 the solution correspond to the extremal solution discussed in section 3. Above 
the shock wave along X~ the solution should agree with the classical shock solution. This 
amounts to setting 

= ^ 2{M -l+x-""- Mx+) ^^-'^ 
as x~ — oo. In practice a large negative initial value of x~ is chosen. 



Appendix B. Numerical method for the RN model 

Here the numerics are also performed in conformal gauge ( [A.l| ) . Asymptotically as 
x~^ — x~ ^ oo the solutions approach the vacuum 

(P = -log{l{x+ -X-)), p = 0. (B.l) 

The equations to be solved are 

d+d-p 



1 _ !lp2(f> 

(B.2) 



d+d-(P = d+d-p + d+(Pd-(P + ^e^p^^'^ 
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Here the boundary conditions are that along = the solution match onto the quantum 
corrected extremal solution [§], while above the shock wave along X~ the solution agree 
with a classical shock solution. This means 

as x~ — oo. 
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